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ABSTRACT 

We present the first survey of submillimeter 12 CO (J = 4 —>■ 3) emission in 
the Magellanic Clouds. The survey is comprised of 15 6' x 6' maps obtained 
using the AST/RO telescope toward the molecular peaks of the Large and Small 
Magellanic Clouds. We have used these data to constrain the physical conditions 
in these objects, in particular their molecular gas density and temperature. We 
find that there are significant amounts of molecular gas associated with most of 
these molecular peaks, and that high molecular gas temperatures are pervasive 
throughout our sample. We discuss whether this may be due to the low metallic- 
ities and the associated dearth of gas coolants in the Clouds, and conclude that 
the present sample is insufficient to assert this effect. 

Subject headings: Magellanic Clouds — galaxies: ISM — submillimeter 
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1. Introduction 

With their proximity, their unextinguished lines of sight, and their profuse star forma- 
tion, the Magellanic Clouds are some of the best extragalactic objects in which to study the 
relationship between molecules and star formation. Because the interstellar medium (ISM) 
in the Magellanic Clouds is deficient in heavy elements and dust, molecular observations of 
these objects probe an interesting regime, perhaps more similar to the conditions in early 
protogalaxies rather than those prevalent today in the Milky Way. Indeed, studies of active, 
metal-poor, nearby dwarf galaxies such as the Magellanic Clouds should offer insight into 
the processes at work in primeval sources. 

What are the effects of low metallicities on the star-forming molecular ISM? We know 
that molecules are more difficult to form both in the gas phase (less O, C, and N) and on 
grain surfaces (fewer grains), and that they are easier to destroy (diminished dust shielding 
of the UV radiation) in low metallicity environments. Thus, molecules other than H 2 are 
more rare in these sources which in particular translates into a dearth of CO emission in 
the Magellanic Clouds and other dwarf irregulars (e.g., Israel et al., 1986; Israel, Tacconi & 
Baas 1995, Taylor, Kobulnicky, & Skillman 1998; Leroy et al. 2005). Furthermore, because 
the FIR and submillimeter lines of the different forms of carbon and oxygen (C + , C, O, and 
CO) dominate the cooling of the star-forming ISM (e.g., Le Bourlot et al. 1993; Wolfire et 
al. 1995), the lower abundances of these elements will make the cooling of the molecular gas 
less efficient. However, at the same time smaller dust-to-gas ratios will yield lower heating 
of the molecular gas, as photoelectric ejection of electrons from small dust grains is the 
chief mode by which starlight heats the gas phase of the ISM. If molecular gas temperatures 
were considerably affected by the metallicity of the ISM, we expect important consequences 
for star formation in such environments. In particular, if the Jeans criterion is relevant to 
star-formation, the mass of collapsable clouds grows for decreasing metallicity as the Jeans 
mass increases Mj ~ T 3 / 2 . Such change could have important effects on the Initial Mass 
Function of stars in these systems. To a first approximation, models suggest that lowering 
the metallicity causes a similar decrement in both heating (by diminishing the dust-to-gas 
ratio) and cooling (by diminishing the C and O abundances; Wolfire et al. 1995). 

It is important to realize, however, that there are many possibilities likely to complicate 
this simple picture. For example, if the dust-to-gas ratio were to decrease faster than 
the metallicity (as suggested by Lisenfeld & Ferrara 1998), or if there were a lack of very 
small dust grains in the low metallicity ISM (as suggested by the faintness of the polycyclic 
aromatic hydrocarbon emission observed toward some of these sources; e.g., Madden 2000), 
the heating processes may become less efficient and the balance may be shifted toward lower 
temperatures. Furthermore, because it is also necessary to consider the metallicity threshold 
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Fig. 1. — 12 CO ( J = 4 — > 3) maps obtained toward our detected sources in the LMC/SMC. 
The coordinates of the central positions of the maps (offsets 0,0) are listed in Table 1. The 
contours are logarithmically spaced, starting at f.5 K km s~ x and increasing by factors of 
1.3 (except for N 79, N 83A, N 113, and N 159W, where they start at 2.5 K km s" 1 ). These 
maps were produced by convolving the data with a Gaussian of FWHM=1'. Most maps are 
sampled on a 30" grid, except for N 48 and N 55A where the grid spacing used was 60". 
The positions marked by the crosses are those of the CO (J = 1 — > 0) sources, where we 
measured the intensities used in the LVG analysis. 
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below which hitherto secondary heat sources (e.g., chemical heating) become important, the 
effects of metallicity on the heating and cooling balance of molecular clouds are very difficult 
to address from a purely theoretical approach. 

Answering some of these questions observationally requires studying the physical con- 
ditions of the molecular gas in nearby, low metallicity sources. At distances of 55 and 63 
kpc, the proximity of the Magellanic Clouds affords single-dish millimeter-wave observations 
excellent spatial resolution attainable in other galaxies only through the use of interferome- 
ters, which permits detailed studies of individual clouds instead of ensemble properties. In 
particular, the ability to spatially separate the emission from different regions makes the 
Magellanic Clouds ideal targets to study the excitation of the molecular gas and its relation- 
ship with star formation. Multitransition studies of CO and other molecules are very useful 
tools to determine the physical conditions of the H 2 (e.g., Johansson et al. 1998; Heikkila 
1998; Chin et al. 1998; Heikkila et al. 1999), but their application is limited if there is no 
information on the intensities of the higher CO transitions, which are extremely sensitive to 
density and temperature. Because of their southern declination, however, there is a dearth 
of submillimeter observations of the Clouds. 

We present here a survey of 12 CO ( J = 4 — > 3) in the molecular peaks of the Magellanic 
Clouds, and the results of an excitation analysis using these and lower J observations. These 
pointings were selected among the brightest 12 CO ( J = 1 — > 0) peaks found by SEST obser- 
vations (Israel et al. 1993), many of which are associated with star-forming complexes and 
Henize (1956) Ha nebulosities, and are thus denoted using the corresponding "N" number. 
In section §2 we present the observations, in §3 we discuss the LVG analysis and its results, 
and in §4 we summarize our conclusions. 

2. Observations and Results 

We observed the (J = 4 — > 3) transition of carbon monoxide ( 12 CO) at v ~ 461.0408 
GHz (650.7 jum) using AST/RO, the Antarctic Submillimeter Telescope and Remote Ob- 
servatory located at the Amundsen-Scott South Pole base (Stark et al. 2001). The ob- 
servations were obtained on the austral winter of 2002, using the lower frequency side of 
the dual AST/RO SIS waveguide receiver (Walker et al. 1992; Honingh et al. 1997), with 
system temperatures of T sys ~ 2000 K. The backend was the 2048 channel low resolution 
(1.07 MHz resolution, 0.68 MHz channels) acousto-optical spectrometer (Schieder, Tolls, 
and Winnewisser 1989). The spectra were observed in position switching mode, chopping 
25' in Azimuth (which is the same as R.A. at the pole). At 461 GHz, the telescope beam 
was measured to have a HPBW~ 109". The forward efficiency determined from skydips was 
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Fig. 2. — Spectra of 12 CO (J = 4 — > 3) emission in LMC sources. These spectra were 
obtained toward the positions marked by crosses in Fig. I. All sources were detected. The 
emission in N 159W is shown scaled down by a factor of two for display purposes. The thick 
black line near the bottom of the plots indicates the velocity and width of the CO emission 
in the lower J transitions, as observed using SEST. 
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~ 70%, and is assumed to be identical to i] mb . Maps of size 6' x 6' centered on the SEST 
coordinates (c.f. Table 1) were obtained for each region, using a 30" grid (in a few cases the 
maps were done on a 60" grid). The data were calibrated using the standard procedure for 
AST/RO, which includes sky, ambient and cold load measurements every 20-30 minutes, 
and processed using the COMB astronomical package. The individual maps are shown in 
Fig. 1, and representative spectra are shown in Figs. 2 and 3. Because the pointing accuracy 
is estimated to be ~ 1', we have selected the emission peaks closest to the center of the map 
within those margins to measure the integrated intensities compiled in Table 1 (indicated 
by the crosses in Fig. 1). Note that this method may introduce a bias in the direction 
of obtaining larger CO (J = 4 — > 3)/(J = 1 — >• 0) ratios. Nonetheless we feel that this 
methodology is justified, as in the Milky Way the positional coincidence between the peaks 
of both transitions is very often observed. In some cases it is apparent that the structure of 
the sources is complex (e.g., N 159W; Bolatto et al. 2000), and in particular for N 167 and 
N 83 the brightest emission peak is well away from the central position. Given the large and 
variable effects of the atmosphere at submillimeter wavelengths and the pointing accuracy of 
the telescope, we estimate the overall absolute calibration accuracy of the data to be ~ 30%. 

Table 1 summarizes our observations. Column 4 lists the 12 CO (J = 1 — > 0) integrated 
intensities observed by the SEST Magellanic Cloud Key Programme toward our sources, 
while column 5 shows the intensities convolved to the angular resolution of AST/RO. Finally, 
column 6 lists the 12 CO (J = 4 — > 3) integrated intensities used in our LVG analysis, 
with their statistical errors. These intensities are those measured in the AST/RO pointings 
marked with crosses in Fig. 1. 



3. Analysis and Discussion 
3.1. Modeling of CO 



The available observed 12 CO and 13 CO line ratios have been modelled using the large- 
velocity gradient (LVG) radiative transfer models described by Jansen (1995) and Jansen 
et al. (1994). These models provide line intensities as a function of three input parame- 
ters: gas kinetic temperature Tk, molecular hydrogen density n(H.2) and CO column density 
per unit velocity (N(CO)/dV). By comparing model line ratios to the observed ratios we 
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Fig. 3. — Spectra of 12 CO (J = 4 — > 3) emission in SMC sources, obtained toward the 
positions marked with crosses in Fig. 1. Only N 12, N 27, and N 83 were detected, and the 
signals are in general considerably weaker than in the LMC. As in the previous figure, the 
CO ( J = 1 — ► 0) velocity and linewidth are indicated by the thick black line at the bottom 
of the plots. 
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Table 1. AST/RO and SEST CO observations of Magellanic cloud H II regions 



Source 


Coordinates 


Observed 12 CO Integrated Intensities in K km s 1 




R.A. 


Dec. 


(43") 


(109") 


(109") 




(B1950) 


(B1950) 


(J = 1 - 0) 


( J = 1 - 0) 


( J = 4 -> 3) 




LMC 


N 48 


05 h 25 m 46!6 


-66°17'36" 


3.7 




6.6±0.7 


N 55A 


05 h 32 m 3O!0 


-66°29'21" 


17.4 


4.9 


2.5±0.5 


N 79 


04 h 52 m 09 ! ?5 


-69 28'21" 


21.3 




6.0±0.7 


N 83A 


04 h 54 m 17 s 


-69°16'23" 


21.1 


6.6 


10.4±1.4 


N 113 


05 h 13 m 40!2 


-69°25'37" 


23.4 




9.7±0.9 


N 159W 


05 h 40 m 01!5 


-69°47'02" 


57.0 


45.6 


20.8±2.2 


N 167 


05 h 44 m 17!6 


-69°23'19" 


19.8 


15.5 


5.6±1.1 


N 214C 


05 h 42 m 21!8 


-71°20 / 33" 


9.2 




2.0±0.4 


LIRL 648 


05 h 14 m 07!0 


-69°38'57" 


13.1 




6.7±1.2 


SMC 


N 12 


00 h 44 m 5O!5 


-73°22'33" 


9.0 


3.46 


3.85±0.7 


N 27 


00 h 46 m 32!9 


-73°21 / 50" 


11.6 


5.63 


5.0±1.3 


SMCB1#1 


00 h 43 m 42!0 


-73°35'10" 


4.1 


2.67 


< 0.5 


SMCB2#6 


00 h 46 m 28 ! ?l 


-73°34 / 10" 


6.4 


3.20 


1.7±0.9 


N 66 


00 h 57 m 26!5 


-72°26'36" 


0.9 


3.41 a 


1.7±1.1 


N 83 


01 h 12 m 29!2 


-73°32 / 40" 


5.6 


2.66 


2.0±0.4 



a Obtained from the 13 CO (J = 1 — > 0) map convolved to a HPBW=109", assuming a 
12 CO/ 13 CO ratio w 11. 

Note. - Nominal map center coordinates are indicated. Column 5 contains the 12 CO 
( J = 1 — ► 0) integrated intensity in the ( J = 4 — > 3) beam for sources with ( J = 1 — ► 0) maps. 
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Table 2. CO line intensity ratios of Magellanic Cloud Objects 

Source 12 CO Transition Ratios 12 CO/ 13 CO Isotopic Ratios References 

(2 -> 1/1^0) (3 2/1 ->•{)) (4 -> 3/1^0) (l->0) (2 1) 



LMC 



N 55A 
N 83A 

N 159W 0.90±0.25 

N167 1.16±0.23 



SMC 



N 12 1.20±0.40 1.0 1.11±0.28 11 9 1, 7 

N 27 0.95±0.35 0.75 0.89±0.30 17 11 1, 8, 9 

SMCB1#1 0.72±0.12 0.4 <0.2 11 13.5 1, 10, 11 

SMCB2#6 ••• ••• 0.53±0.28 ••• ••• 1 

N 66 1.30±0.35 1.0 0.50±0.32 11 7 1, 9, 10 

N 83 1.19±0.31 ••• 0.56±0.22 10 9.5 1,12 



Note. — Ratios are determined from values in Table 1, and from previously published papers, mostly in 
the ESO Key Programme. References: (1) This Paper; (2) Israel ct al. (1993); (3) Israel ct al. (2003); (4) 
Johansson et al. (1998); (5) Bolatto et al. (2000); (6) Garay et al. (2002); (7) Chin ct al. (1998); (8) Heikkila 
et al. (1999); (9) Rubio ct al. (1996); (10) Heikkila (1998); (11) Rubio et al. (1993); (12) Bolatto et al. (2003). 



0.51±0.17 11 ••• 1,2 

1.58±0.52 9 ••• 1,3 

0.9 0.46±0.15 9 5 1, 4, 5 

0.36±0.12 12 ••• 1,6 
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determine the physical parameters best describing the conditions in the observed source. 
In principle, with two isotopes we need to measure five independent line intensities in or- 
der to fully determine the conditions of a single molecular gas component (i.e. 7k, n(i7 2 ), 
iV( 12 CO)/dV, 7V( 13 CO)/dV and a beam filling-factor). By assuming a fixed isotopical abun- 
dance [ 12 CO]/[ 13 CO] = 40 (Johansson et al. 1994) we may decrease this requirement to four 
independent line intensities. As Table 2 shows, this minimum requirement is met by two 
out of nine LMC objects and five out of six SMC objects. The physical conditions of the 
remaining eight objects are, in principle, underdetermined. 

More realistic, and consequently more complex, models of gas excitation that include 
more than one component require many more observations to be properly constrained. Full 
modelling of a two-component molecular cloud using two isotopes requires ten independent 
measurements, which again are reduced to eight by the introduction of a fixed isotopical 
abundance. As this is more than we have actually observed in any of the LMC or SMC 
clouds, it is clear that the solutions may not be unique. As long as the range of possible 
solutions is not excessive, however, they are still useful to constrain the physical parameters 
governing the observed emission. 

We identified acceptable fits by searching a grid of model parameter combinations (10 
K < T k < 150 K, 10 2 cm" 3 < n{ H 2 ) < 10 5 cm~ 3 , 6 x 10 15 cm' 2 ( km s" 1 )" 1 < N (CO) /dV < 
3 x 10 18 cm -2 ( km s -1 ) -1 ) for model line ratios matching the observed values. Although errors 
in the line ratios increase the range of possible solutions, these ratios tend to define reasonably 
well-constrained regions of parameter space. The solutions are somewhat degenerate, as 
variations in the parameters may compensate one another. For instance, a simultaneous 
increase in kinetic temperature and decrease in H 2 densities (or vice versa) yields similar 
line ratios (see §3.2.3). 

3.2. Single— component fits 

3.2.1. Objects with two measured intensities 

For three objects (LMC-N 48, LMC-N 214C, and SMCB2#6) we only have intensities 
in the (J = 1 — >• 0) and (J = 4 — > 3) transitions of 12 CO. The parameters of these clouds 
are thus poorly constrained, and not summarized in a table. Assuming a single molecular 
gas component, we find for SMCB2#6 no effective constraints: T k = 20-150 K, n( H 2 ) = 
10 2 - 10 4 cm" 3 , and N(CO)/dV = 10 16 - 10 18 cm- 2 ( kms" 1 )- 1 . All 12 CO transitions should 
be very optically thick. LMC-N 214C is slightly better determined: temperatures below 30 
K are not allowed and densities appear high n( H 2 ) = 10 4 — 10 5 cm -3 . The ( J = 1 — ► 0) 
transition should not be very optically thick (although the higher transitions are) and its 
isotopic intensity ratio should be 10 - 20, as is in fact commonly observed in the LMC 
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(Israel et al. 2003). Finally, the high 12 CO (J = 4 -> 3)/(J = 1 -> 0) ratio exhibited 
by LMC-N 48 are, in the single-component approximation, only consistent with low optical 
depths (isotopic intensity ratios 20 - 40), rather high densities n( H 2 ) = 10 4 — 10 5 cm -3 and 
temperatures T k = 60-150 K, and gradients N{CO)/dV = 10 16 - 10 17 cm" 2 ( kms" 1 )- 1 . 

3.2.2. Objects with three measured intensities 

For five objects, all in the LMC, (N 55, N 79, N 83 A, N 113 and LIRL 648) we have 
measured 13 CO ( J = 1 — > 0) intensities in addition to the 12 CO ( J = 1 — > 0) and ( J = 4 — >• 3) 
intensities. As is clear from the previous discussion, the molecular gas parameters are still 
underdetermined, but not fully unconstrained (Table 3). Both N 79 and N 113 fit very well 
to a hot and fairly dense model cloud with Tk = 100-150 K, n( H2) = 3000 — 5000 cm~ 3 
and a gradient of about 6 x 10 17 cm~ 2 (kms -1 ) -1 . LIRL 648 is fit, but not very well, by 
a gas at the somewhat lower temperature of 60 K and the somewhat higher density of 
10 4 cm~ 3 . The physical parameters of the molecular gas cloud associated with N 83A are 
inconsistent with the assumption of a single component. Only a very poor fit is obtained in 
the high-temperature, high-density, and high-velocity gradient limit. 



3.2.3. Objects with four or more measured intensities 

Seven objects have a sufficiently large number of measured line intensities to allow a full 
determination of physical parameters, assuming they can be properly described by a single 
molecular gas component. This appears indeed to be the case for SMCBl^l, where we find 
excellent agreement between the observed ratios and those of a model gas characterized by 
T k = 50 K, n(H 2 ) = 700 cm" 3 , N(CO)/dV = 5 x 10 16 cm" 2 ( kms" 1 )- 1 . To fine-tune the 
solutions, as well as to gain insight into the various trade-offs, we have run a finer grid of 
models for this particular source. We summarize the best solutions for the temperature range 
Tk = 30-60 K in Table 4. The overall best fit is achieved using a kinetic temperature Tk = 
50 K, and a density n( H2) = 700 cm~ 3 . At constant temperature, the uncertainty in density 
is about 20%. From the Table, it is also clear that we may allow for a similar temperature 
uncertainty of 20% if we simultaneously increase or decrease the density by 40%. 

The good fit between model and observed line ratios is of particular significance because 
the solution is overdetermined with six independent intensity measurements. SMCBl^l is 
a small and isolated molecular cloud in the SMC Bar, not associated with a star-forming 
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Table 3. Model physical parameters - single component fits 



Source 


Kinetic 


Volume 


Column 




Temperature 


Density 


Density 




T k 


n(H 2 ) 


iV(CO)/dV 




(K) 


( cm~ 3 ) 


( cm _2 km _1 s) 




LMC 


N 55 


30-60 


3xl0 3 


1 x 10 17 


N 79 


100 


5xl0 3 


6-10 x 10 17 


N 113 


150 


3xl0 3 


6 x 10 17 


N 159W 


150 


10 3 


3 x 10 17 


LIRL 648 


60 


10 4 


3 x 10 17 


SMC 


SMCB1#1 


50 


7xl0 2 


0.5 x 10 17 



Table 4. Model fits for SMCB1#1 

Tkin n(H 2 ) N{CO)/dV 7 12 CO// 13 CO 12 CO Rotational Ratio 

(K) (cm- 3 ) (cm- 2 (kms- 1 )- 1 ) ( J = 1 — *■ 0) (J = 2 — > 1) (1-0/2-1) (3-2/2-1) (4-3/2-1) 



Observed constraints 







11.2 


13.5 


1.4 


0.6 


<0.24 


Modeling results 


60 500 


6 x 10 16 


10.4 


12.0 


1.41 


0.56 


0.23 


50 600 


5 x 10 16 


11.4 


12.2 


1.45 


0.54 


0.21 


50 700 


5 x 10 16 


12.0 


12.2 


1.4-3 


0.56 


0.23 


50 800 


5 x 10 16 


12.4 


11.8 


1.40 


0.57 


0.24 


40 1100 


5 x 10 16 


12.4 


11.2 


1.37 


0.58 


0.24 


30 2000 


4 x 10 16 


14.1 


10.8 


1.30 


0.59 


0.24 
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region (Rubio et al. 1993; Reach, et al. 2000; Rubio et al. 2004), thus a low density of ~ 700 
cm -3 is not necessarily surprising, although the kinetic temperature is higher than we would 
have expected for such a cloud in the Milky Way. Table 4, however, illustrates the difficulty 
in pinning down the physical conditions even in a simple cloud with a simple model. 

The measured intensities of N 159-W can also be fit, albeit somewhat poorly, by a single 
hot (Tk = 150 K) and moderately dense component (n( H 2 ) = 1000 cm -3 ). However, none of 
the other five objects have intensities consistent with a single component. Typically, a com- 
ponent fitting the observed 12 CO line ratios would fail completely to explain the 12 CO/ 13 CO 
isotopical ratios, indicating that in addition to a warm gas component the presence of a 
second cooler and dense component should be assumed. The lack of associated luminous 
objects probably explains the exceptionally homogeneous nature of SMCB1#1, required for 
such a well-determined single-component fit to be valid. 

3.3. Dual-component fits 

Since (with the exception of SMCBl^l and perhaps N 159-W) none of the sources 
for which a lot of constraints are available allows a good fit with a single gas component, 
we suspect that the succesful single-component fits of the sources in Table 3 result primarily 
from insufficient information rather than from a simple physical structure. Attempts at more 
sophisticated modelling achieve little of value for most of these sources. The single excep- 
tion is LMC-N 83A where only two observed ratios nevertheless conflict with every single 
component model tried. A dual component model yields acceptable but, not surprisingly, 
poorly constrained solutions (see Table 5). For all the other well-observed sources where 
single-component model fits failed, we have also constructed dual-component fits which are 
likewise listed in Table 5. This table shows that usually one of the two components compo- 
nent is reasonably well-determined, whereas the parameters of the other are generally less 
tightly constrained. 

4. Summary and Conclusions 

A general result of our survey is the detection of significant 12 CO ( J = 4 — > 3) emission 
in most molecular peaks in the Magellanic Clouds. By itself this demonstrates the widespread 
occurrence of significant amounts of warm molecular gas, as tentatively suggested by Israel 
et al. (2003). Application of LVG model calculations show that molecular gas kinetic 
temperatures as high as T kin = 100 - 300 K frequently occur. Detailed analysis of the objects 
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for which multiple line ratios are available strongly suggests that the higher temperatures 
occur in cloud regions that are not very dense (nn 2 = 10 2 — 10 3 cm~ 3 ), and that this gas 
is generally associated with colder (typically T^ iQ = 10-60 K) but much denser (nn 2 — 
10 4 — 10 5 cm -3 ) molecular gas. The simplified analysis possible in cases where fewer line 
ratios are available strongly suggests that the situation in those objects is similar. We note 
that recently reported observations of the N 44 complex in the LMC by Kim et al. (2004) 
show also the presence of strong 12 CO ( J = 4 — > 3) emission, leading the authors to conclude 
the likely presence of very high densities (n Ha ~ 10 5 cm~ 3 ). 

A few clouds are notable in our survey. Already mentioned is the case of SMCBl^l, 
a quiescent, compact molecular cloud devoid of any star forming activity. It appears to be 
essentially homogeneous and of modest density (n H2 = 700 cm -3 ), but is surprisingly warm 
(?kin = 50 K). As there are no known embedded heating sources, this temperature must be 
maintained by the environment of the cloud, located in the south west region of the SMC 
Bar. The other two notable objects also occur in the SMC. N 12 is likewise located in the 
southern part of the Bar. Here, both the dense and relatively tenuous molecular phases 
appear to be surprisingly hot (T^ n = 150 K). In N 66, a very hot (Tkm = 300 K) and a cooler 
(^kin — 40 ±20 K) phase coexist, both at a rather high density of the order of nn 2 — 10 4 cm -3 . 
This particular molecular cloud is a relatively small remnant in a large and luminous star- 
forming complex (c.f., Rubio et al. 1996). It is almost certainly in an advanced stage of 
destructive processing; high densities and temperatures are consistent with such a situation 
(Rubio et al. 2000; Contursi et al. 2000). 

Have we found signs of metallicity effects on the temperature equilibrium of the ISM 
in metal-poor environments? Given the measurement uncertainties and the dearth of com- 
parable datasets on "normal" clouds in our own galaxy such finding cannot be asserted on 
the present measurements. It is suggestive, however, that high temperatures seem pervasive 
even in largely quiescent clouds such as SMCB1#1. This result should be taken with cau- 
tion, as we suffer from a strong sample selection bias: the brightness of the optically thick 
12 CO (J = 1 — > 0) transition is proportional to the source temperature and beam filling 
fraction, and we have pointed toward the brightest 12 CO (J = 1 — ► 0) clouds. Thus it may 
not be surprising that our clouds are warm. Furthermore, studies of Galactic photodisso- 
ciation regions (e.g., S 140, Draine & Bertoldi 1999; NGC 2023, Draine & Bertoldi 2000) 
find that H 2 in star-forming clouds is frequently warmer than the theoretical expectation. 
Similarly, the intensities of the mid-J CO lines tend to be underpredicted by homogeneous 
PDR calculations (Hollenbach & Tielens 1999 and references therein). It is unclear whether 
these discrepancies are due to the presence of other heating mechanisms or to shortcomings 
in the models (e.g., Draine & Bertoldi 1999). 
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From the observational standpoint, a conclusive study of the effects of metallicity on 
gas temperature in the Magellanic Clouds must await the availability of more powerful 
submillimeter instruments observing the southern sky such as the Atacama Large Millimeter 
Array (ALMA), and the recently deployed single-dish telescopes: the Atacama Pathfinder 
Experiment (APEX) and the Atacama Submillimeter Telescope Experiment (ASTE). 

We wish to thank the AST/RO group, and the anonymous referee. ADB wishes to thank 
D. Hollenbach for discussions on the subject of warm molecular gas, and for suggestions that 
helped improve this manuscript. ADB and CLM acknowledge support from NSF grants 
AST-0228963 and OPP-0 126090 respectively. 
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Table 5. Model physical parameters - dual component fits 



Source 


'Cold Dense' Component 


'Hot Tenuous' Component 


Relative 




Kinetic 


Volume 


Column 


Kinetic 


Volume 


Column 


J=2-l 




Temperature 


Density 


Density 


Temperature 


Density 


Density 


12 CO 




Tk 


n(H 2 ) 


N(CO)/dV 


T k 


n(H 2 ) 


N(CO)/dV 


Emission 




(K) 


(cm- 3 ) 


( cm^km" 1 s) 


(K) 


(cm" 3 ) 


( cm- 2 km~ 1 s) 






LMC 


N 83A 


10 


10 5 


10 x 10 17 


100 


10 5 


3 x 10 17 


3:2 




60 


10 5 


1 x 10 17 


150 


5-10xl0 2 


6-10xl0 17 


4:1 




150 


10 5 


3 x 10 17 


150 


10 2 


6-10xl0 17 


2:3 


N 159W 


20 


10 5 


1 x 10 17 


100 


10 2 


1 x 10 17 


1:1 


N 167 


20 


10 4 


0.6 x 10 17 


30-60 


1-I0xl0 2 


0.3-10xl0 17 


2:1 


SMC 


N 12 


150 


10 5 


10 x 10 17 


150 


1-5 xlO 2 


6 x 10 17 


1:1 


N 27 


30 


10 5 


0.6 x 10 17 


60-100 


10 2 


0.6 x 10 17 


1:1 


N 66 


20-60 


10 4 -10 5 


6 x 10 17 


300 


10 4 


0.3 x 10 17 


1:9 


N 83 


10-30 


10 4 


1-2 xlO 17 


100 


3xl0 3 


0.3 x 10 17 


1:8 



